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Abstract 

Wishart  [1961]  and  Takacs  [1963]  derived  the  joint  distribution  of  the  size  N  and  residual 

service  R  as  encountered  by  a  new  arrival  into  a  regular  M/G/l.  Wishart  obtained  the 

OO  •  A  i 

following  expression  for  the  generating  function  Il(z,x)  =  HJ(x)z'7  ,  where  11^  =  Fr[R<x, 
N  =  j) : 

oo 

»<-*>  -  V-*)1-.,)K(1',){  • 

0 

We  exploit  the  fact  that  the  system  size  N  is  known  and  find  the  conditional  r.v.’s  *R;  = 
residual  service  seen  while  N  =  j.*  Our  method  seems  better  suited  to  numerical  work  and  we 
extend  it  to  some  variants  of  M/G/l:  M/G/l/K,  and  then  M/G/l  and  M/G/l/K  with  state- 
dependent  service  and  arrival  rates. 
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Notation  (for  browsers;  the  symbols  are  also  defined  in  context) 

N  =  size  of  system; 

N,  =  size  of  system  provided  server  works; 

A  =  (poissonian)  arrival  frequency 
A  ■  =  arrival  frequency  when  N  = j 
x  =  service  time; 

Xj  =  (state-dependent)  service  time  of  a  customer  whose  service  starts  when  N  =  j  ; 

B(t)  =  Pr(x  <  t); 

®>(t)  =  Pr(x;-  <t); 

P  j  =  Pr(N  =  j); 

P.  =  Pr(server  busy); 

V>(x)  is  an  arbitrary  function  of  x;  ^>(x,y)  is  an  arbitrary  function  of  x,y; 
D2V’(x,y)  =  d\l>(x, y) j dy\  partial  derivative  with  respect  to  the  second  argument 
f^  =  frequency  of  jumps  “i  — ►  j” 

c*x  =  Xi+x2+  +xt  where  the  x,  are  free  copies  of  x  (i.i.d.  r.v.’s  equivalent  to  x 
and  also  independent  of  other  variables  within  the  argument) 

V'(Z)  is  an  arbitrary  function  of  Z  such  that  Eif>( Z)  <  oo  (all  terms  entering  analysis 
are  required  to  have  finite  expectation,  e.g.  E4'\ Z)) 

Erl>( Z)  is  called  the  omni-transform  of  V’(Z).  When  V’(Z)  =  exp(  — sZ)  and  Z>0 

then  Etl’(Z)  is  the  L-S  transform  of  Z 
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Section  1  Service  Residues  as  Conditioned  on  System  Size  in  M/G/l 

Find  R; ,  the  residual  service  time  (“residue”)  of  the  ongoing  service  as  seen  when  N  =  j,  j  >  1. 

Assume  the  P  •  =  Pr(N  =j)  to  be  known  (e.g.  Gross  and  Harris  1985). 

Definition  Let  i /(X)  be  an  arbitrary  function  of  the  process  X  with  finite  Exp(X) 
and  Exl''(X).  We  call  Eip(X)  the  omni-transform  of  X;  if  V’(X)  —  e-s^  we  get  the  L-S 
transform.  The  balance  of  the  process  ip(X)  is  the  equation  EdrlfX)  =  0  for  a  random 
dt. 


The  essence  of  the  omni-method  is  to  study  the  balance  of  t/>(X)  rather  than  of  X.  Among  the 
method’s  advantages  are  freedom  to  choose  i/>,  notational  ease  in  handling  sums  and  mixtures 
of  r.v.’s,  and  bypassing  L-S  transforms  in  many  contexts. 

Definition  Z;  —  Rj  if  N  =  j  and  =  0  if  N  ^  j 

The  process  varies  by  aging  and  jumps:  j—  1  — » j ,  j  +  l-*j,  j  —* j  —  1  and  j-*j-f  1.  We 

assume  Ed^(Z^)  =  0  for  a  random  dt  and  work  out  the  balance.  For  a  random  dt: 

(a)  aging  adds:  FdV-(Z,)|0„ny  =  E[V'(Ri  -  dt)  -  ^(R1)]  =  -dt  P,  E^(RX) ;  dRx  =  -dt 

(b)  jumps  add:  Edt^Z^lo-n  =  dt f01  £|t/-(x)- V'(0)];  f0]  =  AP0 

(c)  jumps  “2-*  1”  add:  JE7d V’(Zx ) 1 2  — ►  i  =  dt f21  £ft(’(x)  —  V’(0)];  f21  =  A2  =  ^Pi 

(d)  jumps  “l-*0”  add:  EdV'(Z1)|1_0  =  dt  f10  £|^(0)  —  tf>(0)]  =  0 

(since  Z,  =  Rj  =  +0  just  before  “1  -»0”  and  Zj  =  0  in  any  state  other  than  “1") 

(e)  jumps  “l-»2”  add:  Edv'^ZJIi  —2  =  dt  f12  E[t/’(0)  —  R-x )] ;  f12  =  APj 
From  (a)  through  (e)  we  get  the  balance  of  V’(Z]),  i.e.  Edv'^Zj)  =  0, 


P1Etl'(R1)  +  U2  g[V( RJ-VKO))  =  (for  -h  £~[V(x)  -  V(0)] 


(1.1) 


where  f01  =  f10  =  AP0  and  f12  =  f21  =  APj .  The  right  side  of  (1.1)  is  known. 


page  3 


Definition  The  omni- convention  calls  for  mentally  applying  the  expectation  operator  E  to  each 
side  of  an  omni-equation;  in  case  of  ambiguity  we  retain  E.  (This  convention  is  kin  to  sum¬ 
mation  convention  in  matrix  and  tensor  calculus.)  E.g.  (1.1)  becomes 


P^Rii  +  f.aMRil-^O)]  =  (foi  +  fai  Mx)-V'(O)] 


(1.2) 


Let  us  consider  the  changes  in  Ei{’(Zj)  for  a  j  >  2.  During  a  random  dt: 

(A)  aging  adds:  Edip(Zj)\asing  =  P ;  E\rJ'{Rj  -  dt)  -  V'(Rj)]  =  -dtP^EV^R,);  dR;  =  -dt 


(B)  jumps  “j-l-*j”  add:  £d^(ZJ-)|i_i_+J-  =  £!v'(R-y_i)  ~  V!(0)];  =  AP;_j 

(C)  jumps  “j  +  l-*j”  add:  £7d^>(Zi)| i+1_ j ;  =  dtfj+1>j  £ty>(x)  -  V’(0)];  fJ  +  1 ;- =  f;  i+1  =  AP_, 

(D)  jumps  uj-*j-l”  add:  £<i V(z>)l>— *>_,  =  E[\f{0)  -  t/>(0)]  =  0 

(E)  jumps  “j-j  +  l"  add:  EdV>(Z^)|  j_J+1  =  dtf;j+1  (V’(0)  —  V’(R-j)]?  fjj+i  =  AP; 

From  (A)  through  (E)  we  get  the  balance  of  t^(Z;) 


p j  tf'(Rj)  +  fjj+i  MRj)  -  V'(Q)]  =  tj-u  (V-(Rj,i)  ~  <KQ))  +  [VW  -  V>(Q)] 


(1.3) 


where  j  =  and  f;  =  f>+i,;  =  A;P^ .  The  right  side  of  (1.3)  is 

known  when  we  solve  for  successive  values  of  j  >  1. 


From  (1.2)  and  (1.3)  we  get  equations  for  moments;  or  L-S  transforms;  or  tail  distributions  of 
Rj  by  setting  ^(R;)  =  R;-‘  for  i  >  1;  or  t/>(Rj)  =  exp(-sR;);  or  t/>(R;)  =  where  =  1  if 
Rj  >  t  and  (j  =  0  if  R;  >  t  for  then  Hy(t)  =  Pr(R;  >  t)  =  E{ . 

Note:  We  can  set  Et/)(R;)  =  Pr(R;  >  t)  if  we  know  that  Pr(R;  >  t)  =  E£(R;)  for  some 
£(R;).  In  linear  omni-equations  with  constant  coefficients,  as  in  our  paper,  we  can  view  V 
as  a  general  functional.  We  need  not  then  find  a  ((R})  and  need  no  omni-convention. 

Let  rp(R})  =  H;(t)  4  Pr(R^  >0)  in  (1.1)  and  (1.2)  and  let  V(xj)  =  Pr(x^  >  t)  =  B(t) .  Then 
,  ,.  ^(R;-dt)-  V'(Rj)  Pr(Ry  — dt  >  t)  —  Pr(R  >  t) 

*  <Ri)  = hm - -  = lim  — - rar - = 
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and  we  get  with  B(t)  =  Pr(x>t) 


-P,  fl{(t)  4-  AP,  H](t)  =  (AP0  +  AP1)B(t) 


(1.4a) 


-P,H^(t)  +  APJHj(t)  =  APJ.1H;.1(t)  +  APJB(t) 


(1.4b) 


H;(t)  can  be  found  from  (1.4a);  from  (1.4.b)  we  can  derive  H  -(t)  if  H;_j(t)  is  known.  Thus  we 
can  find  the  Hj(t)  for  successive  j.  Equations  (1.4)  are  easily  verified  for  M/M/1  with 
Hj(t)  =  B(t)  =  e-p*  for  each  j  .  Moreover  (1.4a)  implies 

j>l  H^(t)-*0  as  t-*oo,  and  P;  H^(t)-»  —  AP;-1  as  t-*0  (1.5) 

From  (1.1a)  and  (1.2a)  we  get  the  recursive  relations  for  R^ 

AP,  Ri  =  pP0  -  (l-p)Pi  and  APJ  +  1  RJ+1  =  AP;R j  -  (1  -  p) PJ  +  J  (1.6) 

If  we  know  all  the  conditional  t£(R;)  we  can  get  the  V’(w).  Clearly 

w0  =  0  and  w;  =  V'(w|N  =  j)  =  V;((j  - 1)  *  x  +  R;),  j  >  1  (1.7) 

V>(w)  =  P0  V’(0)  4-  Pj  V'(Ri)  +  P2  V'(l*x  +  R2)  +  P3  V’(2*x  +  R3)  +  P4  V’(3*x  +  R3)  +  -f  (1.8) 

where  c*x  =  Xj  +  4- xc ;  the  x,  are  free  copies  of  x  (i.i.d.  copies  of  the  generic  x  and 
independent  of  the  Ry). 

Note:  Using  the  renewal  relation  (Krakowski  1987) 

V-(Z)  -  V’(0)  =  2  Et/A'ftZ) ,  Z>0  9k  Z  =  residue  of  Z 
we  can  recycle  (1.2)  and  (1.3)  into  (omni-convention  still  holds!) 


Pt  V’(R:)  4-  fi2  R-i  V'ffi  Rt)  =  (fpi  +  f2) )x  V’(gt.x) 


(1.9a) 


P, gR~T+  f;ij  +  1  Rj  U(^  Rj)  =  fy-i.y  R  j.,  Rj-i)  +  +  |  (1.9b) 

Equations  (1.9)  have  no  derivatives  t/'1  and  thus  in  a  sense  are  integrals  of  (1.2)  an  d  (1.3). 

But  since  both  R;  and  ^R^  arc  arguments  in  (1.9)  there  is  no  labor  saved  unless  we  take 
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a  special  interest  in  the  Sis  R;  in  addition  to  . 


Conjecture :  R;-»  residue  of  xasj-*oo 


page  6 


Section  2  N,  R  in  M/G/l  with  State  Dependent  Service 


We  modify  M/G/l  as  follows.  A  service  which  starts  when  N  =  j  lasts  ,j  >  1.  Our  problem 
is  to  find  the  R;  for  j  >  1  assuming  that  the  P  are  known  (Harris  1967,  Gross  and  Harris 
1985,  pp. 289-292;  Krakowski  July  1986;  a  closely  related  vacation  model  was  treated  by  Harris 
&  Marchal  1988.) 

Let  7jj  ==  Rj  if  N  =  j  and  2,-  =0  if  N  ^  j;  t/-( Z ^ )  is  arbitrary  except  for  Ev'(Z;)<o c  and 
E4<\Zj)  <  oo  .  During  a  random  dt 

(a)  aging  adds:  Edf(Z1)\agtn}  =  Pj  EK'fRj  -  dt)  -  v'-(Ri)]  =  -Pi  dt  Ev'tRj);  dRj  =  -dt 

(b)  jumps  “0-*r’  add:  Edv'(Zi)|0_,  =  dt  f0j  Ej^Xj)  -  v'’(0));  f01  =  AP0 

(c)  jumps  “2-*l"  add:  JE7d V*(Zj )J2 _ ►  t  =  dtf2,  EfV^xJ  —  V(0)];  f21  =  f12  =  APj 

(d)  jumps  “1 -*0”  add:  EdV^Zj)!,  _  0  =  dt  f10  E)V(O)  —  V(0)]  =  0 

(e)  jumps  “1 -»2”  add:  Edt./;(Zi)|,_  2  =  d*  f12  £[^(0)  —  V’(Ha )J ;  fI2  =  AP: 

From  (a)  through  (e)  we  get  (mind  the  omni-convention!) 


Balance  of  *(Zj)  P,  ^(Rt)  +  f„  [l'(Ri)  -  V(0)1  -  (f0i  H-  f21)[^(x,)  -  t/>(0)] 


(2.1) 


where  f01  =  110  =  AP0  and  fJ2  =  f2I  =  Af1,  .  The  right  sid'’  of  (2.1)  is  known. 

For  j  >  2,  during  a  random  dt: 

(A)  aging  adds:  £aV(Z,-)|fl#jn,  =  P;  E[V(R;  -  dt)  -  V’(R,))  =  -dt  P,  V'(R,) !  dR;  =  -dt 

(B)  jumps  “j-l-*j”  add:  £dt/>(ZJ)|J_,_t>  =  dtf^,^  [V( Rj-i )  —  V:(0)];  f,_,w  =  AP;_, 

(C)  jumps  “j  +  l-*j”  add:  Edi/'(Z;)|;  +  1  ,  =  dtfJ+M  E[0(*j)  ~  V’(0)] ;  f;  +  1<;  =  f;,J  +  1  =  AP_, 

(D)  jumps  “j-*j  -  1”  add:  Ed HZj)\ }  _  =  dt  f>,>_ilV’(0)  -  V(0)]  =  0 

(E)  jumps  “j-j  +  l”  add:  Edv»(Z>)l>  —  >+i  =  dt  fJiJ+1  E[V’(0)  -  V’(R;-)) ;  f;,J+1  =  AP; 

From  (A)  through  (E)  we  get  the  balance  of  V'( Z  ■)  (mind  the  omni-convention!) 
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Pj  V^(Rj)  4-  fiJ+1  MRj)  -  V'(0)]  =  fj-i.iWRj-,)  -  V’(0)1  +  fj+1.,  lvA(x,)  -  V-(0)] 


(2.2) 


where  f jj_1  =  fy_j  y  =  Ay_jyPy_j  and  fy  J  +  1  =  fJ+1  y  =  A;P;.  The  right  side  of  (2.2)  is  known 
when  we  solve  for  successive  values  of  j>l. 


When  all  service  lengths  Xy  are  free  copies  of  a  generic  x,  i.e.  V’(xj)  =  V’(x)  for  each  j  >  1,  then 
(2.1)  and  (2.2)  become  (1.2)  and  (1.3)  respectively  -  sis  they  should. 

We  get  Hy(t)  =  Pr(R;  >  t)  by  setting  £V’(R;)  =  Pr(Ry  >  t). 


Note:  It  is  enough  to  know  that  ££(Ry)  =  Pr(Ry>t)  for  some  function  £(R;)  -  we  do  not 
have  to  actually  determine  £((Rj). 


It  follows  that 

j  =  l 


P1H,1(t)-AP1H,(t)  =  AP0-A(P0  +  P,)D1(t) 


j>2 


P;H'(t)-APJHJ(t)  -  AP  _ j 


—  APy_,  Hy_j(t)  —  AP^  By(t) 


(2.3a) 

(2.3b) 


The  right  side  of  (2.3a)  is  known;  so  is  the  right  side  of  (2.3b)  for  each  j  when  the  Hj(t)  are 
known  for  i<j.  (2.3a,b)  imply  that  for  each  j>  1, 


p;  h;(o)  =  a  p 


(2.4) 


From  (2.1)  and  (2.2)  we  get 

APifi-i  =  (Po  +  Pi)Ax,  -  Pj  and  AP;  +  !  R  j  +  ;  =  AP;R;  -  (1  -  Ax;  +  1 )  P;  +  1 
When,  for  each  j  >  1  ,  B;(t)  =  B(t),  then  (2.3a,b)  become  (1.4a,b),  as  they  should. 


j=l 

j>2 


-Pi  Ai(t )  +  AP]  11,(1)  =  ( AP0  +  APj) B(t) 


-P;li;(l)  +  AP;H;(t)  =  AP j_j  H;_,(t)  +  APJB(t) 


(2.5) 

(1.4a) 

(1.4b) 


The  question  arises.  Can  we  derive  the  load  or  delay  from  (2.1)  and  (2.2)?  Unfortunately,  we 
see  no  fair  way.  In  our  model  with  state-dependent  service  the  virtual  load  and  dela>  (and 
kindred  time  lengths)  depend  on  future  arrivals;  this  makes  them  essentially  more  complex. 
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Section  3  The  Rj  and  the  Load  In  M/G/l/K 

Consider  an  M/G/l/K,  i.e.  where  N  <  K;  customers  arriving  while  N  =  K  are  lost.  YVe 
consider  the  load  (backlog,  unfinished  work)  L;  L  =  w,  the  virtual  delay,  when  N<K  but  for 
N  =  K  w  is  not  defined.  We  assume  the  P  -  to  be  known  (Gross  and  Harris,  p.  279-285,  1985). 
Cleariy 

V-(L)  =  P0V-(0)  +  P1^(R,)  +  P2V(1*x  +  R2)  +  P3tf(2*x  +  R3)  +  +  P  K  -  l)*x  +  R3)  (3.1) 

Dexine  now  V’(Z;)  as  before:  «/’(Zj)  =  i/'(Ft;)  if  N  =  j  and  V'(Z;)  =  i/'(0)  otherwise.  The  balance 
of  </'( Z;)  clearly  yields  the  same  equations  for  j  =  1  and  for  1  <  j  <  K  as  the  balance  for  regular 
M/G/l.  Forj  =  K  during  a  random  dt 

(a)  aging  adds: 

£7dv'( ZK)\agin}  =  PK£lv(RA  -dt)-^(RA)]  =  -dtPA.^(RA.);  dRA  =  -dt 

(b)  the  jumps  “K  —  1  -»K”  add: 

Edrli ZA)|A_._  A  =  ^ K 


(c)  the  jumps  “K-*K  — 1”  add: 

Ed^( ZA)|  A_  K_1  =  dtfA_  A_;  £lv'(0)  -  V(0)]  =  0 

Therefore 

j=l:  P.^RJ  +  fnftHRi)-^)]  =  (f0i  +  f2i)[V’(x)-  tf’(O)]  (3.2a) 

1  <  j  <  K:  Pj  ip'(Rj)  +  fJi>+]  [V(R y )  -  v'-(0)]  = 

=  f>-i.i[v,’(Rj-i)-V-(0)J  +  fj+1.J-(V(x)-V'(0))  (3.2b) 

J  =  K  PAv'(RA)  =  fA_jAI^RA_i)-V’(0)]  (3.2c) 

The  tail  distributions  H;(t)  ==  Pr(R;  >  t)  satisfy 

j=l  -PJH;(t)-AP,fl1(t)  =  (AP0  +  AP0)B(t)  (3.3a) 

l<j<K  -Pj»>-i(t)  +  AP,  ftJ-(t)  =  APJ_,HJ_1(t)  +  AP;B(t)  (3.3b) 
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j  =  K  PK  H'A-(t)  =  APk_i  HA._i(t)  (3.3c) 

Starting  with  j  =  l  we  can  compute  the  successive  H^(t).  Defining  R  =  0  for  N  =  0  (readers 
may  favor  other  definitions)  we  have 

V’(N)  R)  =  P0V’(0,0)  +  P1V>(1,Ri)  +  P2V’(2,R2)  +  +  P^K^)  (3.4) 

For  the  load  L  (unfinished  work,  backlog)  we  have 

rp( L)  =  P0  tA(0)  +  Pj  ^(Ri)  +  P2  V’(x  +  R2)  +  P3  ip(2*x  +  R3)  +  +  Pk^((K  — l)*x  +  RK)  (3.5) 

which  appears  to  be  new.  The  economic  motivation  must  be  strong  to  numerically  evaluate 
the  moments  or  the  distribution  of  L;  but  the  complexity  appears  to  inhere  in  the  problem, 
not  in  our  method. 
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Section  4  M/G/l  With  State- Dependent  Service  And  Arrival  Rates 


We  extend  now  M/G/l  with  state-dependent  service  lengths  of  Section  2  by  allowing  also 
state-dependent  arrival  rates.  The  service  time  is  still  determined  at  the  instance  in  which  a 
service  begins  and  equals  x{ ,  i  being  the  size  of  the  system,  including  the  customer  to  be 
served.  The  arrival  rate  is  A;  when  system  size  is  N=j.  Let  us  act  as  if  the  P}  are  known.  So 
soon  someone  derives  these  Pj  our  derivation  of  N,  R  for  this  model  will  be  completed. 

Let,  as  in  Section  2,  Zj  —  R;  if  N  =  j  and  Z j  =  0  if  N  ^j;  >s  arbitrary  but  subject  to 

EipC^j)<oo  and  £'t/’,(ZJ)  <  oo . 


It  is  easy  to  see  that  equations  (2.1a)  and  (2.2a)  stay  valid  but  their  frequency  coefficients 
reflect  the  changed  arrival  and  service  rates.  Clearly  we  now  have 


P1  rp'(R1)  +  f12  MRO  -  HO)}  =  (*oi  +  *21 )Mx,)  -  V'(0)]  (2.1a)=(4.1) 


Pj  4\Rj)  +  fj.j  +  iWRj)  ~  V’(O)]  =  fj-uWRj-i)  ~  m)  +  fj+u  [V’(Xj)  -  4(0)]  (2.2a)=(4.2) 


with  f0J  —  A0  P0 ;  f12  _  f21  —  Aj  Pj ;  f;-_j j  —  A;_j  P;_j ;  fj  2+1  —  fJ+i,j  —  A^P ; . 

From  (4.1)  and  (4.2)  we  get  H;(t)  =  Pr(R;  >t)  by  setting  E4’(Rj)  =  Pr(Rj>t).  Thus 

j  =  l 


P ,  H  j  (t)  —  AP ,  H,(t)  =  A0  P0  —  (A0P  p-l-AjP,)  B](t) 


j>2 


~P;H'(t)-f  AjPjH/t)  =  A;_1P;-1H;-i(t)-fA;P;B;(t) 


(4.3a) 

(4.3b) 


The  right  side  is  known  for  (4.3a),  and  for  (4.3b)  if  we  solve  for  successive  j>l.  (4.3a)  and 
(4.3b)  imply  that  for  each  j>l 


H'(0)  =  0 


(4.4) 


(4.1)  and  (4.2)  imply 


Ai  Ri  —  (Aq  Pq  +  *>i)xi  —  Pi  Qn<*  Pj  Rj  —  ^;-l  **j-l,>  R;'-l  +  A;  x;  P;  (4-5) 
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Section  5  M/G/l/K  With  State- Dependent  Service  and  Arrival  Rates 


We  extend  M/G/l/K  with  state-dependent  services  of  Section  4  by  allowing  state-dependent 
arrival  rates;  service  time  x,-  still  depends  on  system  size  (including  customer  served)  when 
service  starts.  The  arrival  rate  is  A;  when  N  =j.  We  act  as  if  the  P;  are  known,  but  there 
seems  to  be  no  analytical  derivation  yet  available. 

Define,  as  in  Section  2,  Z;  =  Rj  if  N  =  j,  and  Zj  =  0  if  N^j;  j>l.  Let  \p{Z-)  satisfy 
E4'(Zj)<oo,  Erl>\Z-)<  oo,  and  be  otherwise  arbitrary;  then 

V"(Z j )  =  (1  —  P j)^(0)  +  P;  ^(Rj)  (5.1) 

It  is  clear  that  the  equations  (3.3a,b,c)  are  still  valid  so  that 


j  =  l 

Pi  V-'(Ri)  +  f,2WRi)  -  V>(0)]  =  (f01  +  fa,)Mxi)  -  tfr(O)] 

(3.3a)=(5.2a) 

1  <  j  <  K 

Pj^(Rj)  +  f>ti+1[V'(Ri)-V-(0)]  = 

(3.3b)=(5.2b) 

=  f i-u  m] + M*i)  - 

j  —  K 

PA‘  ^(P-a  )  =  fA-7,A  MRA-l)  ~  ^(°)) 

(3.3c)=(5.2c) 

o** 

II 

O 

II 

po  !  Az  =  *21  =  ^l  pi  Pj— i,i  =  fj,j-i  =  ^;-i  P jf— l  ?  Pj.y-n  =  fj  +  i,>  = 

W 

With  the  modified  fi;- ,  the  tail  distributions  H;(t)  =  Pr(Rj  >  t)  satisfy 

j  =  l 

-PjftJCO  +  A^jH.Ct)  =  (AoPo  +  A^jS^t) 

(5.3a) 

1  <  j  <  K 

—  P j  Pj-i(t)  +  A;  P j  H j(t)  =  A;_j  P H ;_j(t)  -f  P;  B ^(t) 

(5.3b) 

j  =  K 

-PKH'A.(t)  =  AK_1PA-_1HK_1(t) 

(5.3c) 

We  can  compute  all  Hj(t)  forj>l.  Defining  (N, R)  =  (0,0)  for  N  =  0 ,  yields 

V’(N,R)  =  P0V’(0,0)  +  P1V>(1,R1)  +  P2v'-(2,R2)+  +  P*  v'(K,Rk)  (4.4) 

We  cannot  from  the  foregoing  find  (certainly  not  easily)  the  load  or  delay  or  any  such  variable 
because  at  any  instant  these  durations  depend  also  on  future  arrivals. 
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Appendix  Joint  Treatment  of  System  Size  and  Residual  Service  In  M/G/l 


We  derive  now  a  single  global  omni-equation  equivalent  to  equations  (1.2)  and  (1.3). 

Definition  V"N,Z)  =  V(0,x)  when  N=0  and  ^N,Z)  =  \£>(N.,R)  when  N  >  1 ;  V>(N,Z)  is  defined 
at  any  time  point. 

Let  N.  =  system  size  provided  N>1.  R  is  defined  only  when  N>1  so  it  needs  no  asterisk.  We 
assume  that  ^(N,Z)  is  a  balanced  r.v.,  so  that  Edrl<(N,  Z)  =  0.  Of  course  we  have 

lKN,Z)  =  P0V’(0,x)  +  P.V>(N.,R)  P0  =  Ax  P.  =  1-Ax  (A.l) 

Let  us  consider  the  balance  of  ^(N.Z)  during  a  random  dt: 

(a)  aging,  which  goes  on  only  while  N  >  1,  adds 
£d^N,Z)|a,inj=P.£;[V’(N.,R— dt)-V’(N.,R)]  =  -dt  P.  ED2  V-(N. ,  R) ;  P.  =  Pr(N>l) 

(b)  arrivals  “0-*l”  add  Ed^N^)^!  =  f01  E(V>(l,x)  -  V-(0,x)J;  f01  =  AP0 

(c)  arrivals  while  N  >  1  add  Ed^(N,Z)|a,  =  f.  £f^(Na.  +  l,Ra)  -  t£(Na.  ,R„)];  f.  =  AP. 

A  subscript  a*  says  that  the  r.v.  is  found  by  an  arrival  into  a  busy  system. 

Since  true  poissonian  arrivals  see,  stochastically,  what  a  continuous  or  random  (poissonian) 
observer  sees  (cf.  Wolff  1982),  we  have  t/’(Na.,Ra)  =  t^(N.,R)  and 

f.  £[V’(Na.  +  l,Ra)  —  V>(N0.  ,Ra)]  =  f.  £JV’(N.  +  1,R)  —  t/’(N. ,  R)] ;  f.  =  AP. 

(d)  departures  (perforce  from  a  busy  system)  add 

Edil’{ N,Z)|d  =  AE[V’(Nd,x)  —  ^(N^  +1 ,0)];  departure  rate  =  arrival  rate  =  A 

Subscripts  d  say  that  a  r.v.  is  seen  by  a  just  departed  customer.  (System  size  is  Nd  +  1  just 
before  a  departure,  and  is  Nd  just  after.)  Since  xp( Nd)  =  V’(Na)  =  V’(N),  and  since  Nd  and  x 
are  independent, 

Ed^(N,Z)|d  =  A£[V’(Nd  ,x)  —  V’(Nd  + 1 ,0))  =  AE[V-(N  ,x)  -  V(N  +  1 ,0)] 
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From  (a)  through  (d)  we  get  (mind  the  omni-convention!) 


-P.  Dz  V-(N.,R)  +  f01  [V(l,x)  -  ^(0,x)]  +  f.  [ifr(N.  +  1,R)  -  V’fN.  ,R))  + 

+  A[N,x)-V(N +1,0)]  =  0  (A.2) 

where  P,  =  1  —  P0  =  p  =  Ax  ;  f01  =  AP0  =  A(l  — p);  and  f,=  arrival  rate  while  server  works 
=  AP.. 

With  V’(N)  =  P0  V’(O)  +  P.  V’CN.)  we  get  from  (A.2)  the  equation 
P.D2V>(N.,R)-AP.(V>tN.  +  l,R)-^(N.,R)]  + 

+  AP.[V-(N.  +  l,0)-V>(N.,x)]  =  APo[V-(l,x)-V>(l,0)]  (A.4) 

Equation  (A.4)  does  not  depend  on  the  definition  Z  =  x  when  N  =  0. 
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